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TRACK  INITIATION  IN  A 
DENSE  DETECTION  ENVIRONMENT 


INTRODUCTION 

Track -while-wan  systems  were  first  proposed  for  surveillance  radars  during  the  1950’s. 

If  the  probability  of  detection  per  scan  is  high,  if  accurate  measurements  are  made,  if  the 
target  density  is  low,  and  if  there  are  few  false  detections,  the  design  of  the  correlation  logic 
and  tracking  filter  is  straightforward.  However,  in  a realistic  radar  environment  these  assump- 
tions arc  never  valid,  and  the  design  problem  is  difficult.  This  paper  will  consider  the 
problem  of  track  initiation  in  a dense  detection  environment. 

In  Fig.  1,  there  are  three  tracks,  and  each  track  is  detected  five  times.  While  it  is  ob- 
vious that  there  are  three  tracks  present,  many  tracking  systems  would  initiate  incorrect 
tracks  because  they  only  associate  the  nearest  detection  with  the  predicted  position  of  a 
tentative  track.  Moreover,  the  situation  in  Fig.  1 rarely  occurs;  the  situation  in  Fig.  2 is 
more  common.  Figure  2 is  of  the  same  three  tracks,  however,  several  detections  have  been 
merged  (i.e.,  individual  targets  are  not  resolved),  three  detections  are  missing,  and  two  faise 
alarms  have  been  introduced. 


RANGf 

Fij;.  1 — History  of  five  scans  of  three  tracks 
showing  all  detections  present 


The  optimal  solution  of  such  problems  lias  been  generated  under  ideal  conditions. 
Specifically,  the  maximum  likelihood  solution  has  been  developed  under  the  assumptions 
that  the  probability  of  detection,  the  probability  of  false  alarm,  the  probability  of  target 
resolution  as  a function  of  target  separation,  and  measurement  error  characteristics  are  all 
known  a priori  and  that  all  targets  are  moving  in  straight  lines.  (A  somewhat  similar 
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Fig.  2 — History  of  five  scans  in  which  detections  weie  missed, 
detections  were  merged,  and  false  alarms  occurred 

approach  was  used  by  Stein  and  Blackman  [1];  however,  they  did  not  consider  resolution 
problems.)  Even  if  all  the  above  assumptions  were  true,  the  maximum  likf-lihood  method 
cannot  be  implemented  in  the  foreseeable  future  because  of  the  enormous  computational 
load.  However,  it  is  still  useful  because  it  provides  a standard  with  which  to  compare  al- 
gorithms that  can  be  more  readily  implemented. 

There  are  two  basic  problems  with  the  maximum  likelihood  method.  The  first  prob- 
lem is  fundamental  and  concerns  the  tendency  of  the  maximum  likelihood  method  to  in- 
dicate *wo  tracks  (with  many  unresolved  detections)  when  a single  detection  is  close  to  a 
single  track.  This  problem  was  solved  by  penalizing  tracks  with  unresolved  targets  or  miss- 
ing detections.  A detailed  description  of  the  maximum  likelihood  method  is  given  in  the 
next  section. 

The  second  problem  is  the  computational  load.  Since  search  techniques  cannot  be 
used  to  maximize  the  likelihood  functions  because  of  the  large  number  of  local  maxima, 
the  concept  of  a “feasible  track”  was  introduced,  a feasible  track  consisted  of  a specified 
number  of  detections  lying  within  a specified  distance  of  a straight  line.  Then  the  maxi- 
mum likelihood  of  occurrence  of  each  combination  of  the  feasible  tracks  was  evaluated. 

If  there  were  iV  feasible  i racks  and  one  is  interested  in  up  to  M track  combinations. 

Zfi , ) likelihood  functions  would  need  to  be  evaluated.  For  instance,  if  .V  = 30  and 

M ~ 4,  the  number  of  likelihoods  calculated  is  31930.  This  problem  is  discussed  in  the 
section  entitled  “Calculation  cf  Maximum  Likelihood.” 

A brief  description  of  how  the  radar  data  are  generated  is  given  in  the  section  entitled 
“Parameters  for  Data  Generation  and  Operation  of  the  Simulation,” and  the  results  of  the 
maximum  likelihood  method  applied  to  various  target  geometries  containing  one  to  four 
tracks  and  several  false  alarms  are  given  under  “Results.”  The  final  section, “Conclusions.” 
summarizes  the  results  and  suggests  a practical  solution  that  is  presently  being  investigated. 
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MAXIMUM  LIKELIHOOD  METHOD 

The  maximum  likelihood  method  involves  calculating  the  total  probability  that  a given 
set  of  detections  correctly  represents  a specified  set  of  tracks.  The  probability  of  detection, 
the  probability  of  false  alarm,  the  probability  of  target  resolution  as  a function  of  taiget  sep 
aration,  and  the  measurement  error  characteristics  are  all  taken  into  account  in  the  calcula- 
tion of  the  likelihood.  To  facilitate  the  mathematical  description  of  the  likelihood  method 
the  following  terms  or  definitions  are  used: 

N( > - the  number  of  scans 

Nt  - the  number  of  tracks 

N D = the  total  number  of  detections 

NpA  = the  number  of  false  alarms 

Ny  = the  number  of  missed  detections  associated  with  the  Nr  tracks 

Nq R = the  number  of  detections  involved  in  resolution  problems  (i.e.,  number  of  de- 
tections used  in  at  least  two  tracks) 

* the  total  number  of  tracks  using  the  k- th  detection  which  is  used  in  at  least 
two  tracks 


xit  = the  range  of  the  detection  associated  with  i-th  track  on  the  /-th  scan.  If  there 
is  no  detection  associated  (i.e.,  track  has  a miss  associated),  xt]  = 0 

y,j  = the  predicted  range  of  the  i-th  track  on  /-th  scan-assuming  a straight  line  tra- 
jectory 


All  of  the  above  variables  are  not  independent.  The  following  relationship  holds: 


where 


nd  =nsnt-  nm  - ntr  + ndr+  nfa  , 


NTRih). 


(1) 


The  difference  between  predicted  and  measured  position  is  assumed  to  be  Gaussian 
distributed  with  zero  mean  and  a variance  of  o2.  Thus, 

1 - (x ,,—y ij)2 /2o2 

(2tto2)1/2  c 


p(xij-yij ) = 


(2) 
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For  later  use,  it  will  he  convenient  to  introduce  the  expression 


f(x,i-y,j)  = 


1,  xtJ  - 0 or  x,j  = xl}  i l 


-(xn-y..)2/2o2  lU 
v ‘J  J,i'  • otherwise. 


Assuming  that  the  probability  of  detection  PD  is  known,  the  probability  of  obtaining  the 
specified  number  of  detections  is 


NsNr 


nsnt~nm  nm 
(pd)  (W>d) 


The  probability  of  not  resolving  any  NTR(h)  = N)t  tracks  which  use  a common  detection 
xtt  is  calculated  by  first  ordering  the  predicted  positions,  so  that 

y.  < y-.  y-.  , 

'1  '2  ‘Nk 

where  for  notational  convenience  the  subscript  denoting  the  scan  has  been  dropped.  When 
one  lets  Dl  = y(.  - y(.  , the  probability  of  not  resolving  any  Nk  tracks  is  given  by 

Nh 

PR(x,.)  = n P(Dt)  (5) 

£=  2 

In  this  case  the  probability  of  not  resolving  two  tracks  separated  by  distance  D (discussed 
in  [2] ) is  given  by 

{1  D < 1.7 R 

(2.6 R-D),‘(.9R)  1.1R  <D  < 2.6#,  (6) 

0 D>  2.6/.’ 

wheie  R is  the  3-dB  pulse  width  (range  ceil  dimension).  Furthermore,  the  position  of  xk 
is  the  sum  of  two  random  variables:  one  uniformly  distributed  between  y,  and  ytNf  , and 
the  other  Gaussian  distributed  with  mean  zero  and  variance  a2.  In  the  appendix  it  “is  shown 
that  the  likelihood  of  xk  can  be  approximated  hy 


whe/e 


Ps(xk ) = (e'f  fe/20  2 ) / Max(y,-W  - v,y  yjlita 2), 


ek  = Max  iO,xk-yiN  yii-xK) 


NRL  REPORT  8238 


;s  'bo  distance  from  xk  to  the  nearest  detection  if  xi:  lies  outside  of  the  interval  defined 
Py  the  predicted  positions;  otherwise  is  zero.  Finally,  the  number  and  position  of  false 
alarms  in  the  range  interval  of  interest  RI  is  given  by  the  Poisson  density 


PFa  ( nfa ) 


N 


(M/) 


FA 


— \ Rr 

e 1 


N FA 

(NFAy.  (Rj) 


Ufa  -xr, 
X e 


WfaV- 


(9) 


A'  fa 

where  X is  the  false  alarm  density  per  unit  length  and  the  (R/)  factor  in  the  denomina- 
tor was  due  to  the  fact  that  the  false  alarms  are  uniformly  distributed  in  range. 

In  terms  of  the  previous  expressions  the  likelihood  of  an  NT  track  combination  is  given 
by  the  following: 


L(Nt)  - Pfa(nFa)  ‘ 


fNsNr 


M 


(PD) 


NsNt-Nm 


N 


(1  ~PD) 


M 


(27 TO2)1/2 


nsnt~nm~ntr  nt  ns 

n n f(xij-yti) 

1=1  j-1 


fl  PR(xk)  pe(xk)- 
k=l 


(10) 


The  first  line  represents  the  false  alarm  probability,  the  second  line  represents  the  de- 
tection probability,  the  third  gives  the  measurement  error  probability,  and  the  last  givps 
the  resolution  probability.  The  maximum  likelihood  method  involves  assigning  to  each 
i-th  track  (yielding  predicted  positions  ytJ)  a sequence  of  detections  (x,y ) that  maximize 
the  values  calculated  by  formula  (10). 


As  presently  formulated,  the  maximum  likelihood  method  would  have  trouble  with  the 
target  geometry  shown  in  Fig.  3.  Let  the  n-tuple  (/ 1 , I2,  ..,  /„)  represent  a track,  where  L 
is  the  detection  associated  with  7-th  scan  of  the  track.  In  Fig.  3,  there  are  two  tracks,  (1,1, 
1,1,1)  and  (2, 2,2, 2 ^W);  M represents  a miss,  and  detection  3 on  scan  1 is  a false  alarm.  How- 
ever, the  maximum  likelihood  method  defined  by  (10)  will  yield  the  solution  involving  the 
three  tracks  (1,1, 1,1,1),  (2,2,2,2/if),  and  (3,2,1  MM)-  The  latter  case  is  more  iikely  (as  de- 
Imed  by  (10))  for  the  following  reasons:  the  false  alarm  likelihood  has  been  increased  by 
6 X 104  (by  eliminating  the  false  alarm),  the  detection  likelihood  has  been  decreased  by  0.6 
(12  out  of  15  detections  instead  of  9 out  of  10  detections),  the  measurement  likelihood  is 
increased  by  removing  a (l/27ro2)!  factor  (two  detections  declared  resolutions  but  one  de- 
tection added)  and  by  eliminating  two  Gaussian  errors,  and  the  resolution  likelihood  is 
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Fig.  3 - Scan  history  showing  a 
false  track  created  as  the  result  of 
a false  alarm  and  the  incorporation 
of  detections  of  actual  tracks 


decreased  by  1/27 ro2.  Thus,  since  6 X 104  is  greater  than  (l/27ra2)1 /2,  the  likelihood 
for  three  tracks  is  larger  than  the  likelihood  for  two  tracks.  As  formulated  by  (10),  the 
maximum  likelihood  method  will  always  try  to  eliminate  false  alarms  by  introducing  false 
(extraneous)  tracks. 

To  eliminate  this  problem,  two  factors  have  been  introduced.  One  factor  penalizes 
tracks  that  have  unresolved  detections,  and  the  second  factor  penalizes  tracks  that  have 
missing  detections.  Thus,  the  modified  likelihood  is  given  by 

NTR-NDli 

Lm(Nt)  = L(Nt)(Fr)  (Fm)  . (11) 

Usually,  we  take  1 >FM>FR.  The  values  presently  being  used  are  Fm  =0  .2  and  FR  = 0.1 
For  the  rest  of  this  paper,  the  maximum  likelihood  method  will  be  implemented  by  (11). 
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CALCULATION  OF  MAXIMUM  LIKELIHOOD 

Search  techniques  cannot  be  used  to  maximize  the  likelihood  function  (11)  because  of 
the  large  number  of  local  maxima.  To  solve  this  computational  problem  the  concept  of  a 
feasible  track  was  introduced.  Then  the  maximum  likelihood  of  each  combination  of  feas- 
ible tracks  was  evaluated. 

In  this  study,  five  scans  were  considered,  and  a feasible  track  required  at  least  three 
detections.  Furthermore,  all  detections  in  a feasible  track  were  required  to  be  within  2.6 
range  cells  of  the  line  joining  the  first  and  last  detections.  If  two  feasible  tracks  differ  onlv 
by  misses,  for  instance,  (/]  ,/2,  / 3,  J4,  /5)  and  (A/,.  /2,  /3,  14,  A/),  the  track  with  additional 
misses  is  retained  only  if  its  velocity  differs  from  the  other  track  velocity  by  more  than 
9.14  m/s  (30  ft/s),  where  velocity  is  determined  from  the  first  and  last  detections. 

Next  the  maximum  likelihood  is  calculated  for  all  single  tracks.  A direct  search  tech- 
nique is  used  to  determine  the  target’s  position  on  the  first  and  last  scan.  For  each  detec- 
tion associated  with  the  track  it  is  determined  whether  it  is  more  advantageous  to  label  the 
detection  as  coming  from  the  track  (with  its  associated  Gaussian  error)  or  whether  to  de- 
clare a false  alarm  and  a missed  detection*.  At  the  end  of  this  process,  if  a track  has  a de- 
tection for  each  scan,  it  is  called  a “perfect”  track. 


Next,  the  likelihood  is  calculated  for  each  possible  two-track  combination.  That  is, 
if  there  are  30  feasible  tracks,  there  are  30(29)/2  = 435  two-track  combinations.  However, 
if  the  two  tracks  do  not  have  any  common  detections,  the  two  tracks  are  said  to  be“iso- 
lated”and  the  maximum  likelihoods  for  the  single  tracks  are  used.  If  the  two  tracks  do 
have  common  detections,  each  track  is  considered  to  extend  from  its  first  to  its  last  detec- 
tion. For  each  detection  associated  with  only  one  track,  it  is  determined  whether  it  is 
more  advantageous  to  label  the  detection  as  coining  from  the  track  (with  its  associated 
Gaussian  error)  or  whether  to  declare  a false  alarm  and  a missed  detection.  For  detections 
common  to  both  tracks,  one  of  three  actions  is  determined  upon:  1)  to  declare  the  tracks 
unresolved,  2)  to  associate  the  detection  with  the  nearest  track  and  declare  one  target  miss, 
or  3)  to  declare  two  missed  detections  and  a false  alarm.  It  should  be  noted  that  all  of  the 
previous  calculations  proceed  on  a scan-to-scan  basis.  Therefore,  it  is  possible  to  obtain  a 
slightly  different  likelihood  if  the  scans  were  evaluated  in  a different  order  (e.g.,  if  one 
introduced  a miss  on  scan  2,  one  may  not  want  — or  be  allowed  — to  introduce  a miss  on 
scan  3). 

After  all  two-track  combinations  are  evaluated,  ali  three-track,  four-track,  etc.,  com- 
binations are  evaluated  Usually,  if  the  true  answer  is  an  n-track  combination,  all  n+l-traek 
combinations  are  evaluated.  Next,  the  best  track  combinations  (usually  the  best  5 to  10 
are  saved)  are  maximized  by  the  use  of  direct  search  techniques  in  which  each  target’s  posi- 
tion on  the  first  and  last  scan  is  varied.  Finally,  the  track  combination  with  the  maximum 
likelihood  is  chosen  as  the  correct  series  of  tracks. 


* A miss  is  never  introduced  if  the  miss  lowers  the  number  of  detections  below  'hat  required  for  a 
feasible  track. 
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When  a large  number  of  tracks  is  present,  the  computational  time  on  NRL’s  ASC  com- 
putor  can  become  exorbitant.  For  instance,  calculation  of  the  likelihoods  for  a four-track 
combination  of  .50  feasible  tracks  requires  over  a minute.  Thus,  to  increase  computation 
speed  the  method  was  modified  to  take  advantage  of  “perfect”  tracks  — those  that  have  de- 
tections on  each  scan.  Since  it  is  very  likely  that  the  perfect  track  will  be  in  all  the  high 
likelihood  track  combinations,  only  track  combinations  that  include  the  perfect  track  (or 
tracks)  wili  be  evaluated.  For  instance,  if  there  are  30  feasible  tracks  and  tracks  2 and  8 
are  perfect  tracks,  only  one  two-track  combination,  28  three-track  combinations,  and  (28) 

( 27 )/2  four-track  combinations  will  be  evaluated.  Thus,  for  this  example  only  407  (1  + 

28  + 28( 27 )/2 ) track  combinations  need  be  evaluated  instead  of  all  the  31  900  possible  track 
combinations:  (30(29/2)  + 30(29i  (28)/6  + 30  (29)  (28)  (27)/24). 


PARAMETERS  FOR  DATA  GENERATION  AND  OPERATION  OF  THE  SIMULATION 

Before  the  results  of  several  simulations  are  given  the  data  generation  technique  will 
lie  described  briefly.  The  targets  are  assumed  to  be  travelling  in  straight  lines  at  constant 
speeds.  The  radar  detections  are  generated  on  a scan  basis  in  the  following  manner:  A 
decision  is  made  on  whether  or  not  a target  is  '.elected.  If  a target  is  detected,  its  position 
is  calculated  according  to  the  straight  hue,  and  a Gaussian  wander  is  added  to  its  position. 
Next,  false  alarms  are  generated  according  to  a PoLson  density,  and  all  the  detections  are 
ordered  in  range.  The  detections  are  examined,  and  it  is  decided  whether  adjacent  detec- 
tions should  be  merged.  If  a detection  is  not  merged,  a Gaussian  measurement  error  is 
added.  If  several  detections  are  merged,  all  merged  detections  are  replaced  by  a single  de- 
tection whose  range  is  a Gaussian  measurement  error  added  to  a detection  unirormly  dis- 
tributed between  the  nearest  and  farthest  merged  detections. 


Data  Input  Cards 

The  data  generation  and  simulation  operation  are  controlled  by  four  or  six  input  cards. 
The  parameters  covered  by  the  cards  and  their  formats  are  as  follows. 

Card  1 (15  parameters,  1515  format) 

1.  The  number  of  targets  ( N ) 

2.  The  number  of  scans 

3.  The  first  repetition  case 

4.  The  last  repetition  case 

5.  The  number  of  best  (high-likeiihood)  track  combinations  saved 

6.  The  smallest  track  combination  (minimum  of  1) 
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7.  The  largest  track  combination  (maximum  of  5) 

8.  Starter  for  random  number  generator 

9.  Track  indicator:  either  tracks  inputted  (NCAL  = 1)  or  tracks  calculated  (NCAL.  Nfi.  1) 

10.  Print  control:  IPR  = 2 (mimimum),  1PR  = 1 (intermediate),  or  IPR  = 0 (detailed) 

11.  Number  of  misses  allowed  in  a feasible  track 

12.  If  NREV  = 0 or  1,  do  only  1 power  iteration;  otherwise  do  2 . 

13.  NOPT  = 0 or  optimal  number  of  tracks  (optimal  track  groups  must  be  supp'ied  as  the 
last  cards) 

14.  Is  a minimum  number  of  perfect  tracks  required?  If  negative,  do  not  limit  to  perfect 
tracks. 

15.  IPLOT  = 0,  no  plot;  IPLOT  = 1,  plot  ranges;  and  IPLOT  = 2,  plot  normalized  ranges 

Card  2 (6  parameters, 6K10.2  format) 

1.  Probability  of  detection,  actual 

2.  Average  number  of  false  alarms  in  range  interval,  actual 

3.  Variation  of  true  position  (in  range  cells),  actual 

4.  Variation  of  measured  position  (in  range  cells),  actual 

5.  Average  track  velocity  (ft/s),  actual 

6.  Standard  deviation  of  track  velocities  (ft/s),  actual 

Card  3 (7  parameters,  7F10.2  format) 

1.  Far  range  of  interval  of  interest  (ft) 

2.  Near  range  of  interval  of  interest  (ft) 

3.  Range  cell  dimension  (ft) 

4.  Allowable  velocity  change  of  each  track  from  average  (0  to  1) 

5.  FSVEL.  (IF.  0.1:  no  velocity  constraint  for  individual  feasible  tracks 
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6.  Average  track  velocity  (ft/s),  assumed 

7.  Standard  deviation  of  track  velocities  (ft/s),  assumed 

Curd  d (9  parameters,  9F5.2  format) 

1.  Probability  of  detection,  assumed 

2.  Average  number  of  false  alarms  in  range  interval,  assumed 

3.  Scan  time  (s) 

4.  Variation  of  true  position  (in  range  cells),  assumed 

5.  Variation  of  measured  position  (in  range  cells),  assumed 

6.  Allowable  distance  between  detection  and  track  (in  range  cells) 

7.  Allowable  velocity  difference  between  feasible  tracks  (ft/s) 

8.  Penalty  factor  for  resolutions  U'K ) 

9.  Penalty  factor  for  missed  detections  ( ) 

If  tracks  are  specified  (NCAL  = 1,  parameter  9 on  card  1),  as  opposed  to  being  ran- 
dom, cards  5 and  6 are  used. 

Card  5 ( N parameters,  8F10.2  format) 

1.  Range  (ft)  of  first  target 

2.  Range  (ft)  of  second  target 

N.  Range  (ft)  of  AMh  target 

Card  6 ( N parameters,  8F10.2  format) 

1.  Speed  (ft/s)  of  first  target 

2.  Speed  (ft/s)  of  second  target 

N.  Speed  (ft/s)  of  N- th  target 
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If  one  desires  to  evaluate  a specific  track  combination,  one  can  accomplish  this  by 
setting  NOPT  A 0 (parameter  13,  card  1)  and  by  supplying  one  card  for  each  repetition. 

ImsI  Cards,  one  for  each  repetition  (NOPT  parameters,  1515  format) 

1.  Feasible  track  1 

2.  Feasible  track  2 


NOPT.  Feasible  track  NOPT 


RESULTS 

Information  in  Tables 

The  pertinent  parameters,  excluding  the  target  parameters,  are  given  in  Table  1 , and 
the  target  parameters  are  given  in  Table  2.  The  maximum  likelihood  approach  was  applied 
to  10  independent  realizations  of  the  5 cases  given  in  Table  2;  the  results  are  summarized 
in  Table  3.  Of  the  50  cases  run,  7 were  incorrectly  identified.  However,  it  was  judged 
(by  the  authors)  that  all  7 incorrect  solutions  were  the  most  reasonable  result.  In  most 
cases  where  the  number  of  tracks  was  underestimated,  die  true  track  contained  fewer  than 
three  detections  and  these  were  judged  to  be  false  alarms.  The  two  cases  where  the 
track  had  a velocity  error  greater  than  10%  occurred  because  the  track  either  used  a false 
alarm  or  stole  a defection  from  another  track. 

Examples  — Maximum  Likelihood  Method 

Ten  examples  of  radar  detections  of  a four-track  situation  will  be  reviewed  to  illustrate 
the  maximum  likelihood  method.  The  detections  made  on  five  scans  on  each  repetition  are 
shown  in  Figs.  4(a)  through  4(j).  In  each  figure  the  total  range  interval  is  3.7  km  or  2 n.mi. 
Note  that  for  presentation  purposes  the  ranges  have  been  normalized  by  adding  1524  m 
(5000  ft)  per  scan,  which  corresponds  to  a velocity  of  305  m/s  or  1000  ft/s.  In  the  figures, 
dots  represent  detections,  M’s  represent  missed  detections,  arcs  represent  unresolved  detec- 
tions, and  FA  indicates  a false  alarm.  The  dashed  lines  represent  the  true  tracks.  In  each 
scan  the  detections  are  numbered  from  right  to  left. 
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Table  1 — Simulation  Parameters 


Parameter 

Value 

Number  of  scans 

5 

Number  of  misses  allowed  in  track 

2 

Probability  of  detection 

0.85 

Average  number  of  false  alarms  per  scan 

0.3 

Gaussian  wander,  standard  deviation 

30.5  m (100  ft) 

Gaussian  measurement  error,  standard  deviation 

30  5 m (100  ft) 

Range  interval 

2.0  (n.mi.) 

Range  cell  dimension 

152.4  m (500  ft) 

Scan  time 

5.0  s 

Penalty  factor  for  resolution 

0.1 

Penalty  factor  misses  (FM) 

0.2 

I 


i 


i 

i 


1 


Table  2 — Target  Parameters 


Case 

No. 

No.  of 
Targets 

Initial  Ranges  (km) 

Velocities  (m/s) 

1R 

1 

U(181.6,  185.3)* 

G(305,  15)t 

2 

2 

183.8,  182.9 

305,  282 

3R 

3 

U(181.6  to  185.3)* 

G(305,  15)f 

3 

3 

183.6,  183.3,  182.9 

305,  305,  305 

4 

4 

183.8,  182.9,  182.7 

290,  274,  305,  274 

*U(I81.G,  185.3)  indicates  that  the  initial  target  positions  are  uniformly  distri- 
Ijuted  between  the  two  ranges  given. 

t<«(305,  15)  indicates  that  the  velocities  are  Gaussian  distributed;  the  first  repre- 
sents the  mean  value,  and  the  second  gives  the  standard  deviation. 
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Fig.  4 ~ The  diagrams  (a)  through  (j)  making  up  this  figure  present  10  iterations  of  one 
simulated  4-target  raid  and  the  radar  responses  that  occurred.  The  variations  in  results 
are  caused  by  random  false  alarms  that  were  introduced  and  by  noise  ano  clutter.  Dots 
represent  detections,  M's  represent  missed  detections,  arcs  represent  unresolved  detect- 
ions, and  FA  indicates  a false  alarm.  The  dashed  lines  show  the  true  tracks.  In  each 
scan  the  detections  are  numbered  from  right  to  left. 
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in  each  scan  the  detections  ar?  numbered  from  right  to  left 

15 


TRUNK  AND  WILSON 


SCAN 

NO 


? I 

I 


J 


I 

I 


!• 

I 

I 

I 

I 

M 

I 

I 

I 

♦ 


4 

I 


I 

I 


J 


/ 

I 

/ 

/ 


/ 

r 

/ 


/ 

/ 


• 

i 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

i 


T 

i 

i 

i 

I 

T 


e 


i* 

/ 

/ 


4 

/ 

/ 

/ 

/ 

;• 

/ 


/ 


I / 

1 I 
! / 

■ / 

1 / 

* 


SCAN 

NO 


4 


- 3 / km  <2r»  m*#  - 
RANGE  INTfRVAl 


•/ 

/ 

/ 

i 


i 

» 


4 

/ 


• / 

i , 

i / 
i j 

fr* 


3 7km  ( ? r>  i ) 
RANGf  INTf  RVAl 


0) 


0) 


Fi[>.  4 (Continued)  — The  diagrams  (a)  through  (j)  making  up  this  figure  present  10 
iterations  of  one  simulated  4-targct  raid  and  the  radar  responses  that  occurred.  The 
variations  in  results  are  caused  by  random  false  alarms  that  were  introduced  and  by 
noise  and  clutter.  Dots  represent  detections,  Af's  represent  missed  detections,  arcs 
represent  unresolved  detections,  and  FA  indicates  a false  alarm.  The  dashed  lines  show 
the  true  tracks.  In  each  scan  the  detections  are  numbered  from  right  to  left. 
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In  Fig.  4(a),  the  selection  in  accordance  with  the  maximum  likelihood  is  the  follow- 
ing combination  of  four  tracks:  (1,1,  A/,  A/,1),  (2,2,A/,1 ,1 ),  (3,A/,  A7,2,2),  and  (4, 3,2,3, 3). 
The  nearest  false  likelihood  involving  track  (1 ,1 ,1  ,A/,A/)  instead  of  (1 ,1  ,A7,A7,1 ),  was  only 
five  times  smaller. 

In  case  4(b),  the  three  tracks  ( 1 ,1 ,1  ,A/,1 ),  (2, 2, 2, 1,2),  and  (3, A/, 3, 2, 3)  were  selected. 
Although  incorrect,  this  obviously  is  what  a reasonable  person  would  select.  The  closest 
likelihood  to  this  solution  differs  by  a factor  of  1000. 

In  case  4(c),  the  correct  tracks  (1,1, 1,1,1),  (A/, 2, 2, 2, A/),  (2, A/, 3, 3, 2),  and  (3,3,4 ,4,3) 
are  chosen.  All  other  track  combinations  considered  are  simple  variations  of  the  above- 
tracks. 

In  case  4(d),  the  correct  tracks  (1,1, 1,1,1),  (A/, 2, 2, 2, A/),  (2, 3, 3, 3,2),  and  (A/,4,4,4,3) 
are  selected.  The  closest  likelihood,  which  is  not  a trivial  variation,  differs  by  a factor  of 
10  000. 

In  case  4(e),  the  correct  track  combination  (1,1, 1,1,1),  t2,2,2 ,A/,A/),  (2, 3, 3, 2, 2),  and 
(A/,4,4,3,3)  had  the  largest  likelihood.  The  largest  likelihood  of  a three-track  combination, 
ignoring  track  2,  differed  by  1000. 

In  case  4(f),  the  track  combination  (1,1, 1,1,1),  (2, 2, 2, M, A/),  (A/, 3, 3, 2, 2),  and 
(3, 4, 4, 3, A/)  was  selected.  All  other  combinations  considered  were  simple  variations  of 
this  case. 

In  case  4(g),  the  correct  track  combination  (A/,1 , A/,1 ,1 ),  (A/,2,1 , A/,1 ),  (I, A/, 2, 2, 2), 
and  (2,3, 3, 3,3)  was  selected  even  though  there  were  only  five  detections  on  the  first  two 
tracks.  The  second  largest  likelihood,  which  dropped  track  2,  was  2o  times  smaller  than 
the  maximum  likelihood. 

In  case  4(h),  the  four-track  combination  (1,1 ,1,1,1),  (A/,2,2,2,2),  (2, 3, 3, 3, 3),  and 
(3. 4, A/, 4, 4,)  was  selected.  The  second  track  selected  two  lalse  alarms  (detections  2)  on 
scans  4 and  5 instead  of  detections  1.  The  likelihood  of  the  true  track  combination 
differs  from  the  maximum  by  a factor  of  five. 

In  case  4 (i ),  the  correct  tracks  (1,1, 1,1,1),  (2, 2, 2, A/, A/),  (3, 3, 3, 2, 2),  and  (4, M, 4,3, 3) 
were  selected.  The  likelihood  of  the  three-track  combination  ignoring  track  2 is  100  000 
times  smaller. 

In  case  4(j),  the  correct  tracks  (1,1,1,A/,A/),  (2, A/, 2,1 ,1 ),  (A/, 2, 3, A/, 2),  and  (3, 3,4, 2,3) 
are  selected  and  the  closest  three-track  combination  differs  by  a factor  of  10  000. 

In  summary,  in  the  10  repetitions,  two  false  track  combinations  were  selected. 
However,  both  of  these  were  very  reasonable  solutions.  That  is,  with  the  given  detections 
these  are  the  tracks  one  wouh'  expect  any  operator  or  algorithm  to  deduce. 
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Since  the  maximum  likelihood  solution  assumes  that  the  probability  of  detection 
fPjj),  probability  of  false  alarm  (PpA),  and  Gaussian  measurement  error  (<;),  are  all  known 
a priori,  a sensitivity  analysis  of  the  four-track  c rmbination  was  conducted.  In  the  first 
case,  the  probability  of  detection  was  assumed  1.)  be  0.95  instead  of  the  true  value  of  0.85, 
in  the  second  case,,  tiie  average  number  of  false  alarms  per  scan  was  assumed  to  be  0.6 
inste  ad  of  the  true  value  of  0.3;  in  the  third  case,  the  Gau»sian  error  was  assumed  to  be 
61  m (200  ft)  instead  of  the  true  value  of  30.5  m (100  ft);  and  in  the  la->t  case,  all  the 
incorrect  assumptions  were  made.  The  results  are  shown  in  Table  4.  The  three  repetitions 
that  produced  different  results  were  1,  7,  and  10.  In  case  1 (referring  to  Fig.  4(a))  when  a 
larger  Gaussian  error  was  assumed,  track  (1 ,1 , A/, A/,1)  was  replaced  by  track  (1,1 .1,3/, A/). 
This  had  the  effect  of  removing  a false  alarm.  In  case  7,  different  tracks  were  produced 
when  one  assumed  = 0.95  and/or  o = 61  m (200  ft).  The  resulting  three  tracks  (see 
Fig.  4 (g))  are  (M.2,1  ,M,1 ),  ( 1 ,A/,2,2,2),  and  (2, 3, 3, 3, 3);  that  means  that  track  ,1) 

is  no  longer  detected.  In  case  10,  when  all  the  incorrect  assumptions  were  made,  tra.  k 
(2, M, 2, 1,1)  was  dropped,  resulting  in  only  the  three  tracks  (1,1, 1,1,1),  (2,2,3,37,2),  and 
(3, 3, 4,9, 3)  being  detected.  In  general,  the  maximum  livelihood  method  is  rather  insensitive 
to  the  assumed  parameters.  The  parameter  that  it  is  most  sensitive  to  is  the  Gaussian  error. 


Table  4— Number  of  Tracks  Estimated  for  4— Track  Case 
When  Incorrect  Parameters  Are  Used 
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CONCLUSIONS 

'Hie  maximum  likelihood  method  of  initiating  tracks  works  extremely  well.  However, 
the  method  cannot  he  implemented  because  of  the  enormous  computational  requirement. 
For  instance,  it  took  40  seconds  on  the  NRI,  ASC  computer  to  evaluate  all  possible  four- 
track  combinations  of  45  feasible  tracks.  Thus,  a more  practical  procedure  must  be  con- 
sidered. Presently,  a modification  of  a raid  detector  studied  by  Flad  [3|  is  being  pursued. 
The  basic  idea  is  to  declare  a target  raid  and  estimate  the  raid  velocity  and  number  of 
targets  in  the  raid. 
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Appendix 

LIKELIHOOD  OF  MERGED  TARGETS 


If  several  targets  are  merged,  the  position  of  the  unresolved  detection  is  given  by 

x = y -f  2,  (1) 

where  y is  uniformly  distributed  between  plus  and  minus  A (the  nearest  and  farthest 
predicted  target  positions  of  the  merged  targets)  and  z is  a Gaussian  measurement  error 
with  mean  0 and  variance  rr2.  The  density  of  x is  given  by  the  convolution, 


p(x)  = 


L 


_1 1 

2A  (27 to2)* 


c-{x-y)2l 2ci“'  dy 


(2) 


Equation  (2)  will  now  be  evaluated  for  the  two  special  cases  when  A » a and  a » A. 
If  A » a and  |x|  < A (i.e.,  detection  is  between  predicted  positions),  the  integral  of  the 
Gaussian  density  is  approximately  1,  and  (2)  reduces  to 


p(x)  = 


1 

2 A ' 


(3) 


If  A » o but  |jc|  > A (i.e.,  detection  outside  predicted  positions),  p(x)  is  approximately 
given  by 


where  |xj  = A + 6 
and 


f4) 


<HT)  = 


— e"“2  du. 


(5) 


It  should  be  noted  that  the  situation  |*|  > A will  rarely  occur  when  A » ci. 

When  c » A and  |x|  < A (which  will  be  very  rare),  the  exponential  is  essentially 
one,  and  p(x)  reduces  to 
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P (X)  = 


I 

(27T02)1/2 


(6) 


On  the  other  hand,  when  o » A and  |x|  > A,  the  exponential  is  essentially  constant  and 
can  be  pulled  outside  the  integral,  resulting  in 


P (*)  = 


27T02 


1/2 


e-S  2/2 o 2 


(7) 


where  [x|  = A + 8.  Combining  (3)  through  (7),p  ( x ) can  be  approximated  by 


where 


P (x)  = 


e-i  2/2 a 2 

Max{2A , (27T02)1  /2}> 


(8) 


5 = Max{0,  x-A,  . 
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